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Microscopic models for heavy-fermion materials often assume a local, i.e., momentum- 
independent, hybridization between the conduction band and the local-moment / electrons. Moti- 
vated by recent experiments, we consider situations where this neglect of momentum dependence is 
inappropriate, namely when the hybridization function has nodes in momentum space. We explore 
the thermodynamic and optical properties of the highly anisotropic heavy Fermi liquid, result- 
ing from Kondo screening in a higher angular-momentum channel. The dichotomy in momentum 
space has interesting consequences: While e.g. the low-temperature specific heat is dominated by 
heavy quasiparticles, the electrical conductivity at intermediate temperatures is carried by unhy- 
bridized light electrons. We then discuss aspects of the competition between Kondo effect and 
ordering phenomena induced by inter-moment exchange: We propose that the strong momentum- 
space anisotropy plays a vital role in selecting competing phases. Explicit results are obtained 
for the interplay of unconventional hybridization with unconventional, magnetically mediated, su- 
perconductivity, utilizing variants of large-iV mean-field theory. We make connections to recent 
experiments on CeCoIns and other heavy-fermion materials. 



PACS numbers: 75.20.Hr,75.30.Mb,71.10.Li 



I. INTRODUCTION 

Recent years have seen a revival of research in heavy- 
fermion metals, due to the wealth of fascinating phe- 
nomena which can be found in these materialsjii 2 -^^ 
These include non-trivial charge and spin order, un- 
conventional superconductivity, non-Fermi-liquid behav- 
ior, as well as quantum criticality beyond the Landau- 
Ginzburg- Wilson paradigm. The standard microscopic 
description of heavy-fermion systems is based on versions 
of the Anderson or Kondo lattice models, consisting of 
conduction (c) electrons and local moments on a regu- 
lar lattice, with a spatially local hybridization or Kondo 
coupling between c-electrons and local moments. The 
formation of a heavy Fermi liquid in such a model origi- 
nates from the screening of the local moments at low tem- 
peratures through a lattice generalization of the Kondo 
effect - this is reasonably well understood, e.g., using 
slave-particle or dynamical mean-field approaches. 

For some materials, recent experiments^ 7 - indicate that 
the assumption of a local (i.e. momentum-independent) 
hybridization is insufficient for a full understanding of the 
data. For example, optical conductivity measurements 7 
in CeMIns (M = Co,Ir,Rh) do not show the conventional, 
frequently observed hybridization gap^^ ' 10 ^ 1 ^ 2 but in- 
stead have been interpreted in terms of a distribution of 
gap values. Microscopically, a momentum dependence in 
the hybridization is not surprising, as the local-moment 
orbitals are usually of / type and may hybridize with sev- 
eral conduction-electron orbitals. While in certain cases 
this momentum dependence does not lead to significant 
changes in observable properties (as compared to a local 
hybridization), the physics can be qualitatively different 
if the hybridization has zeros (i.e. nodes) in momentum 
space. In analogy with unconventional superconductors, 



having a pair wavefunction with non-zero internal angu- 
lar momentum, we may call heavy-fermion materials "un- 
conventional" Fermi liquids if the Kondo electron-hole 
pairs have non-zero angular momentum. Such uncon- 
ventional heavy Fermi liquids, formed below the coher- 
ence temperature in the described setting, will have e.g. 
quasiparticles with strong momentum-space anisotropy, 
to be discussed in more detail in this paper. (It is worth 
pointing out that the existence of hybridization nodes 
does not imply that parts of the local moments remain 
unscreened.) 

Importantly, unconventional hybridization will influ- 
ence the entire complex phase diagram of heavy-fermion 
compounds, where the lattice Kondo effect and various 
types of long-range order compete for the same elec- 
trons. Clearly, a strong momentum dependence of the 
hybridzation may favor or disfavor certain ordering phe- 
nomena. For instance, a dichotomy in momentum space 
arising from anisotropic Kondo physics may determine 
which (unconventional) superconducting phase is realized 
at lowest temperatures. 

On the theory side, hybridization with higher an- 
gular momentum has been discussed in a few papers 
only. Refs.Hlil studied the case of a half-filled con- 
duction where a Kondo semimetal replaces the conven- 
tional Kondo insulator - this physics is likely relevant to 
CeNiSn and CeRhSbJ^ Recently, Ghaemi and SenthiliS. 
examined aspects on "higher angular-momentum Kondo 
liquids" , starting from a Kondo lattice model with non- 
local Kondo coupling. Some of their results in the Fermi- 
liquid regime are related to ours below, and we shall 
comment on similarities and differences. Let us note 
that differentiation of electronic properties in momen- 
tum space is a common theme in correlated electron sys- 
tems: For instance, in the copper-oxide high-temperature 
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superconductors, quasiparticles properties are known to 
vary strongly along the Fermi surfaced 

The purpose of this paper is a detailed investigation of 
heavy-fermion metals featuring hybridization functions 
with momentum-space nodes. In the first part, we study 
important Fermi-liquid properties including Fermi sur- 
face, effective mass, specific heat, and optical conductiv- 
ity. We also discuss the temperature-dependent electrical 
resistivity. The focus of the second part is on ordering 
phenomena competing with Kondo screening: Here we 
concentrate on magnetically mediated superconductiv- 
ity. Assuming a direct exchange interaction between local 
moments, we discuss mean-field phase diagrams focussing 
on the interplay of hybridization symmetry and pairing 
symmetry. Most of the concrete calculations are done 
using the slave-boson mean-field approximation on two- 
dimensional lattices, but most of our ideas apply more 
generally, including to situations where the hybridiza- 
tion does not have nodes, but otherwise varies strongly 
in momentum space. 

The remainder of the paper is organized as follows. In 
Sec. |TT] we introduce the microscopic model to be em- 
ployed; its mean-field treatment is subject of Sec. IIIII 
We shall point out the relation between Kondo and An- 
derson lattice models with non-local coupling between c 
and /-electrons. Sec. IIVI is devoted to properties in the 
Fermi-liquid regime, comparing local with non-local (un- 
conventional) hybridization. In particular, the optical 
conductivity will be calculated and discussed in relation 
to experiments on the CeMIns compounds. 7 In Sec. [V] 
we shall touch upon the low-temperature physics beyond 
the slave-boson approximation. Sec. IVII discusses quali- 
tative properties of the electrical resistivity, in particular 
the perturbatively accessible temperature regime above 
the single-impurity Kondo temperature Xk. The compe- 
tition of Kondo screening and superconducting pairing is 
subject of Sec. IVIII Two regimes will be distinguished, 
depending on whether the transition temperature T c is 
comparable to or smaller than Fermi-liquid coherence 
temperature T co h. An outlook will conclude the paper. 



II. MODEL 

In this paper, we shall restrict our considerations to 
two-band models of heavy-electron materials, with sim- 
ple tight-binding hopping of electrons. We shall general- 
ize the Anderson and Kondo lattice models to the case 
of a non-local coupling between the conduction (c) and 
local (/) electrons, and discuss the relation between two 
models. 

The Hamiltonian of the Anderson lattice model is given 

by 

Halm = E( £k ~ ^) c L c k<7 + E( e / ~ rifLf** (*) 

kcr kcr 

+ E ^(/L c kCT + 4,/kJ + u E n M n m 

kcr i 



where c' k(7 (/ k(T ) creates a conduction (/) electron with 
momentum k, spin a and energy ek (e/). /- and o 
electrons are hybridized via the momentum-dependent 
hybridization Vk. Finally, U is the on-site Coulomb re- 
pulsion of /-electrons. The chemical potential /i influ- 
ences the band fillings n c and nj of the c and /-electrons, 
respectively. (Note that Halm has the full translational 
invariance of the underlying lattice.) 

While a full microscopic treatment of the / electron 
lattice would require to consider the low-energy Kramers 
doublet state, e.g., of one / electron in a f 1 configuration 
of Ce in the presence of spin-orbit and crystal-field inter- 
actions, we shall proceed with the simplified model |T]). 
Hence, we shall take the momentum dependence of Vk as 
a phenomenological input, noting that it is dictated by 
the lattice structure and the overlap of c and / orbitals. 
(In principle, the momentum dependence of Vk can be 
further rcnormalized by interaction effects, see Sec. IVl 
below.) The directional dependence of Vk as function 
of k/fc may be expanded into spherical harmonics. Con- 
ventionally, one neglects all non-zero angular-momentum 
components and assumes a local hybridization Vk = V. 
This paper is concerned with situations where the zero 
angular-momentum component is small or vanishes, and 
hence the momentum dependence of Vk can no longer be 
ignored (because e.g. Vk displays nodes in momentum 
space). It is convenient to decompose Vk = V/3k, where 
the form factor /?k is dimcnsionless and normalized to, 
e.g., J2k Pk = N w here Af is the number of lattice sites. 
It is obvious that the thermodynamics as well as most 
other observables of the system will depend on |Vk| 2 only, 
exceptions will be noted in the course of the paper. The / 
level is assumed to be non-dispersive; this approximation 
is relaxed in Sec. IVIII 

In the Kondo limit, i.e., V — > oo, U — * oo, ej — > 
— oo with V 2 /ef finite, charge fluctuations are frozen. 
A Schrieffer- Wolff transformation^ which projects out 
empty and doubly occupied states of the / levels, leads 
to a Kondo lattice model: 

Hklm = E ^ c L c k, + E 2Jkk'e- l ( k '- k)R < S. t • s k , k < 

kcr kk'i 

v. / 

Hj 

(2) 

where s kk ' = c L- c W c l<W 2 ! and = «k - M- 

In the Kondo limit, nf is fixed to unity, and Si is the 
local moment at site i formed out of the /-electrons. 
(An additional potential scattering term arising from the 
Schrieffer- Wolff transformation will be neglected.) To 
leading order, the Kondo coupling is 

J kk < = 2F k Vk/ (j^— + — J = Jo/3k/3k< (3) 
\ u + e f ~ e fJ 

where Jo = V 2 ( u * e + ~r£j)- I n re& l space, the Kondo 
interaction Tij takes the form 

T~t.J = 2Jo E Pn-iPm-iSi ' s mn (4) 
imn 
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where P n -% denotes the Fourier transform of /3k de- 
pending on the distance |R„ — Rj|, and s mn — 



y , 



'ma aa' ^na 



/2 is a non-local conduction-electron 



spin density. 

Importantly, for each impurity Sj, Eq. ^ describes 
a single- channel Kondo model, with the exchange "sym- 
metry" determined by the hybridization symmetry of the 
underlying Anderson model. In contrast, Ghaemi and 
Scnthil— start out from a Kondo lattice model where 
each local moment is exchange-coupled to neighboring 
conduction-electron sites as follows: 



U'r 



: E 



T S 



(5) 



This is a multi- channel Kondo model, where the different 
screening channels correspond to different linear combi- 
nations of the conduction electrons at the surrounding 
sites m. In such a model, the screening channels are of 
different strength, and the strongest will dominate the 
low-temperature physics. The authors of Ref. [TH argue 
that, depending on microscopies (i.e. lattice and band 
structure properties), a higher-angular- momentum chan- 
nel (e.g. d-wave) can dominate over the conventional 
symmetric (s-wave) channel. The general relation be- 
tween Anderson and Kondo models with non-local cou- 
pling has been discussed e.g. in Ref. UtI There it was ar- 
gued that, in an Anderson model, the coupling to a corre- 
lated conduction band opens new screening channels, and 
the effective model will be a multi-channel Kondo model, 
because the charge fluctuations in the band (accompany- 
ing the non-local hopping) are suppressed by conduction- 
electron correlations. 

For our purpose, we note that the two Kondo lattice 
models with TLj and Tl'j (with one screening channel 
dominating) become equivalent in the slave-boson mean- 
field (saddle-point) analysis employed below, because the 
different screening channels correspond to different sad- 
dle points. Therefore, our mean-field results derived 
for the single-channel models (|1I2|) can be directly com- 
pared with the ones for the multi-channel Kondo model 
in Ref. [H 



III. MEAN-FIELD APPROXIMATION 

To obtain quantitative results, we shall employ the 
standard slave-boson mean-field approximation for the 
Kondo and Anderson lattice models (the latter with in- 
finite U). In the following, we briefly summarize the cor- 
responding formalism ^ 18 ' 19 Note that in all cases we re- 
strict our attention to states with spatial translational 
invariance. 

We start with the Anderson model. For infinite on-site 
repulsion, the three states of each / orbital can be repre- 
sented by auxiliary fermions fi a and spinless bosons rj, 
such that the physical / electrons fi a — rjfi a , together 



with the constraint 

V f f f. +rV = 1 

/ j J l(J<> 1(7 1 % % 



(6) 



At the saddle point, the slave bosons condense, (r^) = r. 
which implies a rigid hybridization between the c and the 
/ bands. The mean-field Hamiltonian of the Anderson 
lattice model reads 



Halm.mf = £ £kC L c k CT + e/ZL/k 

kcr kcr 

+E^K/Xck CT +cLJw) 

kcr 

- A ( £ flh. + N{r 2 - 1)) - /, (j2 4 

V kcr / V kcr 




where A is the Lagrange multiplier implementing the con- 
straint ([6]) at the mean-field level, and the effect of the 
chemical potential fi on the / electrons has been absorbed 
in A. The three mean-field parameters r, A, and [i are 
obtained from minimizing the free energy, leading to the 
self-consistency equations: 



Y,WL^ + h.c.) = 2AAAr, 



(8a) 



E<£UJ = m-r 2 ), (8b) 

kcr 

E< c LO = Mn c . (8c) 

kcr 

The expectation values can be easily expressed in terms 
of the Green's function of the diagonalized mean-field 
Hamiltonian, for details see App. [AJ 

For the mean-field analysis of the Kondo lattice model 
one represents the local moments Si by auxiliary fermions 
fia, S l = f Xw Jh a aa'Ji<j>i witn tne constraint 



VP 7 =1- 



(9) 



The Kondo interaction takes the form 

Hj = -J Pn-iPrn-Jlc^clrJ^, (10) 



imnaa' 



where additional bilinear terms have been dropped, as 
they can be absorbed in chemical potentials. The Kondo 
interaction term TL j can be decoupled using auxiliary 
fields b, conjugate to (- J J2 na Pn-if^c na ), i.e., 6, re- 
flects the hybridization between the / and c bands at 
site i. At the saddle point, the 6j condense, and trans- 
lational invariance dictates b., = b. The Kondo lattice 
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mean-field Hamiltonian is 
Wklm.mf = ^ e kC ktT C 



k<r 



kcr 



&EM c Uka 

kcr 



-4 

•A) 




(11) 

As before, the three parameters b, Ao, and /i are deter- 
mined by self-consistency equations which now read: 



kcr 

E(/kJ kCT > = ^. 

kcr 

I>Lo = 



2b 



AA-, (12a) 
Jo 



(12b) 
(12c) 



k<T 



These mean-field equations are equivalent to the ones of 
the Anderson lattice model, Eq. ([8]), if the Kondo limit 
is taken there, for details see App. [Bj 

We note that mean-field Hamiltonians (|7|lljl repre- 
sent the N = oo saddle-point solutions of certain SU(7V) 
Anderson and Kondo lattice models. In this mean-field 
picture, both the Anderson and Kondo lattice models 
are mapped onto two-band systems of non-interacting 
fcrmions with a self-consistently determined renormal- 
ized hybridization between the bands. At high temper- 
ature, the slave-boson condensation amplitude vanishes 
(leading to two decoupled bands), whereas the condensa- 
tion amplitude is finite below the single-impurity Kondo 
temperature Tk- In the Kondo lattice model, Tk is given 
by: 



2 



N ^ e k 2T K 

k 



(13) 



with e k = £k~ M- The neglect of fluctuations in the mean- 
field approach causes an artificial finite-temperature 
transition at 7k which can in principle be cured by in- 
cluding the coupling toaU(l) gauge field, see Sec. IVl 

The effective two-band picture is appropriate to de- 
scribe the low-temperature Fermi-liquid regime, i.e., the 
quasiparticle physics for temperatures below the Fermi- 
liquid coherence temperature T co h. An approximation for 
T co h can be extracted from the T — slave-boson solution: 
For the Kondo lattice model one obtains T co h = b 2 /D 
where D is the bandwidth. 19 (Also, Ao oc T co h-) The 
Fermi surface resulting from the slave-boson approxima- 
tion fulfills Luttinger's theorem: The momentum-space 
volume enclosed by the Fermi surface is proportional to 
the total number of electrons ntot = n c +nf where nj = 1 
in the Kondo limit. At small non-zero temperature T, 
the mean-field parameters acquire quadratic T correc- 
tions characteristic of a Fermi liquid, e.g., [b(T)—b ]/bo oc 
-T 2 /Tl where 6 = b(T = 0). 



The simple two-band picture of the slave-boson ap- 
proach has been confirmed, e.g., using the dynamical 
mean-field theory 2 ^ (DMFT) approach to the Anderson 
lattice model, which fully includes local correlations and 
inelastic processes: The DMFT results 2 ^ nicely show the 
formation of a coherent heavy band crossing the Fermi 
level at low temperatures, rather well separated from the 
second band. 



IV. LOW-TEMPERATURE PROPERTIES OF 
THE FERMI-LIQUID STATE 

This section will discuss the properties of the heavy- 
Fermi liquid state in a heavy-fermion system with un- 
conventional hybridization. Effects of inter-moment ex- 
change and possible ordered phases will be ignored, we 
will come back to these issues in Sec. IVIII Quantita- 
tive results will be obtained by solving the slave-boson 
mean-field equations for different model parameters, but 
the qualitative aspects will be of general validity unless 
otherwise noted. 

We restrict ourselves to two-dimensional systems on a 
square lattice. For the c-electrons a tight-binding disper- 
sion will be assumed: 



£k 



-2t (cos fc T + cos k„ 



(14) 



Results will be shown for hybridization functions of the 
form V4; = VbAc with 



/?k = 




s-wave 
cos k y extended s-wave 
cos k y d x 2_ y 2-w&ve 



(15) 



The d a; 2_ y 2-wave and the extended s-wave case corre- 
spond to different linear combinations of hybridization 
between an /-site and its nearest-neighbor c-sites (also 
discussed by Ghaemi and Senthfli^), while the s-wave 
case with local hybridization is shown for comparison. 

In addition, we will also consider a lattice appropriate 
to model the Celn planes of the CeMIns materials (M=Ir, 
Rh or Co) . Those crystallize in a HoCoGas-type tetrago- 
nal structure ; 22 ' 23 where the Ce and the in-plane In ions 
are located on two interpenetrating square lattices. Thus 
we assume a c electron dispersion as in Eq. (|14p , and hy- 
bridization functions of the form 



Ac = 



2 cos % cos ~y extended S- wave 
2 sin sin D xy -wave 



(16) 



Hybridization functions /3 k which formally break inver- 
sion symmetry, i.e. have odd angular momentum I, are 
possible as well. Observables are governed by V fc 2 and dis- 
play inversion symmetry; hence there is little qualitative 
difference between even and odd angular momentum hy- 
bridization (apart from the location of the hybridization 
nodes). Exceptions are transport anisotropies for / = 1, 
briefly discussed in Sec. IIVEI 
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along the antinodal direction. Note that the function 
fcp (</>), where is the angle parameterizing the k space 
direction, can be multivalued due to the momentum de- 
pendence of Vk: Analyzing the equation Z2k = one 
finds that this generically happens for small band filling 
n c (Fig. [2ji) . In the case of extended s-wave hybridiza- 
tion, Fig. [2b, two Fermi sheets emerge, as the lower band 
crosses the Fermi level twice. 

For the hybridizations of Eq. (| 16|) , arising from two in- 
terpenetrating square lattices of c and / electrons, sample 
Fermi surfaces are shown in Fig. [3] In the d-wave case, 
Fig. [5Jd, the nodal lines are simply rotated by 45 degrees 
w.r.t. to Fig.Hh,d. 



FIG. 1: (Color online) Dispersion relation of the two ef- 
fective bands, 21,2k (|17[) . and the corresponding momentum 
distribution function n^, for a d x 2_ y 2-w&ve hybridization on 
a 2d square lattice, a) along the momentum-space diagonal, 
b) along a direction which encloses the angle 0.1957T with 
the k x axis. The parameters of the Kondo lattice model are 
Jo/t = 2.0 and n c = 0.4. 



A. Band structure and Fermi surface 

To simplify the discussion, we shall work in the Kondo 
limit. The eigenvalues of the mean-field Hamiltonian 
(fTTj) . representing the effective bands of the heavy Fermi 
liquid, are given by 



21, 2k 



-A + e k ± \/(Ao + £k 



(17) 



This band structure is illustrated in Fig. [T] Along cer- 
tain "nodal" lines in momentum space the hybridization 
vanishes, and the two bare bands of the c and / parti- 
cles cross (Fig. QJ,); for d x i_ y i symmetry this applies to 
k x = ±k y . Otherwise, the hybridization causes a band 
repulsion (Fig. [TJa) , which is maximum in the "antinodal" 
direction. As the bare / band is non-dispersive, the two 
bands zi^k do not overlap, and consequently only one 
band crosses the Fermi level. For less than half filling, 
n c < 1, the Fermi surface is thus determined by Z2k = 0. 

Fig. Q] also shows the momentum distribution func- 
tion of the c-electrons, = (cj c c k ). It shows a jump at 
the Fermi wave vector kp , the jump height given by the 
quasiparticle weight Z , sec below. As typical for heavy- 
fermion systems, also shows a rounded step at the 
"small" Fermi surface of the original c electrons, i.e., the 
the Fermi surface in the absence of a hybridization (or at 
temperatures T 3> T co h). 

Sample results from a full numerical solution of the 
mean-field equations (fT2")) are shown in Figs. [5J[3]and in 
Figs. [3J [5] below. Fig. [^displays 2d Fermi surfaces for the 
hybridization functions of Eq. (|15[) , with parameters cho- 
sen such that the specific heat coefficient is identical in 
the four cases, see below. For d x 2_ y 2-wave hybridization, 
Figs. [2t,d, one clearly observes a "small" Fermi momen- 
tum kp along the nodal directions, while kp is "large" 



B. Thermodynamic properties 

The leading low-temperature thermodynamics can be 
directly obtained from the effective two-band description 
of the slave-boson approximation^ For non-interacting 
fermions, the Sommerfeld coefficient, 7 = Cy/T, of the 
specific heat is related to the density of states (per spin) 
at the Fermi level, A = N(w = 0), through 7 = ^-N . 
In standard isotropic Fermi-liquid theory, the effective 
mass m* is defined through the slope of the dispersion 
at the Fermi level, Vp — kp/m*, where vp and kp are 
Fermi velocity and Fermi momentum, respectively. The 
density of states is Ao = m*k F ~ 2 /d in d dimensions 
where C2 — 2tt and C3 = 2ir 2 . 

In anisotropic systems, a suitable definition for a 
(direction-dependent!) effective mass is: 



1 



m*(k) 



<9e k 



(18) 



FS 



where e k is the quasiparticle energy of the band crossing 
the Fermi surface (FS), and k = |k|. Then, the density 
of states is given by 



A 



FS 



27T 



(19) 



Note that Ref. [l5| defined a quantity m* via the second 
(instead of the first) derivative of quasiparticle energies, 
which in general plays only a subleading role in thermo- 
dynamics. 

In our two-band system, ek = 22k for n c < 1. In d = 2 
some factors of A; drop out, such that 



N () = 



1 



2tt 



dA v 



4tt 2 J |V k 2;2k| 

FS 



(20) 



where m* is the effective mass as extracted from specific- 
heat measurements (i.e. the effective mass of an isotropic 
Fermi liquid with the same 7), and dA^ is the Fermi sur- 
face element. Under the condition that the Fermi surface 
can be parameterized in the way k = kp(4>) this integral 
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a) 



(jr,7r) (— g,u) 



{km) (-k,x) 




b) 



c) 





(2r,jr) {— k,k) 
~1 d) 



(jt,.t) 




(7r,-7r) f— 2r,— ?r) 



(Jr.— jf) [— M— m) 



(jt,-jt) i—g—g) 



(jt,-jt) 



FIG. 2: Fermi surfaces for a) s-wave, b) extended s-wave, and c,d) d x 2_ y 2-wa,ve hybridization. The band filling is n c = 0.3 
in panels a-c, whereas n c — 0.1 in panel d; the Kondo coupling is chosen such that the specific heat coefficient is identical in 
all four cases: a) J /t = 1.0 (Ao = -0.025, fx = -2.47, b = 0.281), b) Jo/t = 0.89 (Ao = -0.0067, fx = -2.3, b = 0.281), c) 
Jo/t = 0.97 (A = -0.0061, fi = -2.46, b = 0.173), d) J /t = 1.957 (A = -0.0385, fj, = -3.698, b = 0.466). 



can be rewritten as 



An 
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V(M0)) 2 + (fcM0)) s 

|Vk£2k|fc=fc p (^) 



(21) 



Sample results for the quasiparticle velocity, 
|Vk^2k|FS: as function of the Fermi-surface angle 
(i.e. the direction) are shown in Fig. [3] for the vari- 
ous hybridization cases introduced in Eq. (|15p . The 
microscopic parameters are chosen such that all four 
cases lead to the same value of the density of states 
An and thus the same specific-heat coefficient. The 
corresponding total effective mass is around 125 times 
the bare c-electron mass. For d x 2_ y 2-wave hybridization 
the velocity and therefore also the inverse effective mass 
has a maximum at the nodal line (0 = 7t/4) - here m*(k) 
corresponds to approximately the bare c-electron mass. 
Away from the nodal line the velocity rapidly decreases. 
In contrast, for both s-wave-like hybridizations, the 
velocity is approximately constant (and small) along the 
Fermi surface. 

The electronic quasiparticle weight, Z(k), can be easily 
extracted as well. Z measures the overlap between the 
physical c electron and the low-energy quasiparticle at 



(~K,-n) 



a) 




{-n,-n) 



{n—n) {—n—n) 



{ti—ti) 



FIG. 3: Fermi surfaces using the hybridization functions (|16[) 
for interpenetrating c and / square lattices, a) Extended S- 
wave with J /t = 0.46 (A /* = -0.011, fx/t = -2.38, b/t = 
0.187), b) A^-wave with J /t = 1.975) (A /t = -0.001, 
fi/t = —2.34, b/t — 0.138). In both cases, n c = 0.3, and 
the parameters are chosen such that the 7 coefficient is the 
same as for the data in Fig. [5] 



the Fermi surface. In the mean-field approach of two 
hybridized bands, Z is given by»^ 



Z(k) 



Oik - £k) 2 



(22) 



Results for Z are displayed in Fig.0 It is unity along the 
nodal lines of the hybridization, but becomes very small 
away from it - the latter is the typical heavy-fermion 
situation. For the s-wave-like hybridizations, Z turns out 
to be independent of the momentum direction; in the s- 
wave case this follows from Vk = const, whereas in the 
extended s-wave case this follows from the coincidence of 
the momentum dependence of the hybridization Vk and 
the c electron dispersion ek- 

At this point, let us comment on a few important is- 
sues. First, even in the presence of hybridization nodes, 
all local moments of the Anderson or Kondo lattice are 
fully screened in the low-temperature limit. This is obvi- 
ous from the slave-boson solution which clearly describes 
a Fermi liquid, but also beyond slave bosons we see no 
reason for a (partial) breakdown of Kondo screening: For 
instance, in DMFT complete screening will occur once 
the effective bath density of states at the Fermi level is 
finite (which is the case here). Physically, the local mo- 
ments are entities in real space, whereas the hybridization 
nodes are defined in momentum space. Second, as the 
nodes cover only a set of momenta of zero measure, hy- 
bridization nodes do not easily lead to so-called two-fluid 
behavior (i.e. a heavy-Fermi liquid coexisting with local 
moments), which has been advocated on phcnomcnolog- 
ical grounds^ We note that these statements also hold 
if both quasiparticle bands (c-like and /-like) cross the 
Fermi level. Although one may speculate about the exis- 
tence of gapless spinons at the / Fermi points where the 
hybridization vanishes^ these would again cover only a 
set of momenta of zero measure. 



C. Influence of a magnetic field 

Let us briefly discuss the effects of a weak external 
field applied to the heavy Fermi liquid. In general, a 
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FIG. 4: The quasiparticle velocity at the Fermi level, |VkZ2k|, in a logarithmic plot vs. momentum-space angle cj> for a) s-wave, 
b) extended s-wave, and c,d) d 2 ,2_ H 2-wave hybridization. Parameters are as in Fig. [2^-d. In all four cases, the "total" effective 
mass (as derived from the specific heat) is around 125 times the bare electron mass. 




FIG. 5: As in Fig. [4] but now showing the quasiparticle weight Z in a logarithmic plot vs. angle 4>. 



Zeeman field will cause a spin splitting of the Fermi sur- 
face, with spin- and field-dependent effective masses [or 
densities of states N a (B)]2£- The qualitative field depen- 
dence, N rJ {B) — N (l + aB/B ) with B proportional to 
the Kondo temperature, is not changed by a momentum 
dependence of Vk. (The field dependence of the mean- 
field parameters leads to a subleading correction oc B 2 
to N a (B)££-) The anisotropy of Vk of course causes an 
anisotropy of the fc-space distance between the spin-split 
Fermi sheets, as the Fermi velocity is highly anisotropic. 

For magneto-oscillation measurements, the cyclotron 
mass is an important quantity, given by 27rm c = 
dA(E) /dE, where A(E) is the area enclosed by the quasi- 
particle iso-energy curve in a momentum-space plane per- 
pendicular to the applied orbital field. Thus, in dimen- 
sions d > 2 the cyclotron mass m c depends on the field 
direction. However, in 2d this dependence is absent, and 
m c is identical to the (averaged) quasiparticle mass m* 
extracted from the density of states or specific heat, in- 
dependent of momentum-space anisotropics. 

D. Optical conductivity 

The optical response of heavy-fermion metals has been 
studied extensively.® Experiments probing the optical 
conductivity a(uj) usually show a Drude peak well sep- 
arated from mid-infrared excitations. These features 
have been interpreted in the two-band picture advocated 
above: While intra-band particle-hole excitations pro- 
duce conventional metallic Drude-like response, inter- 
band excitations lead to finite weight at elevated ener- 
gies. In a picture of free fermions, the threshold energy 
of these optical inter-band excitations measures the min- 
imum gap between occupied and unoccupied states in 
the lower and upper band, respectively. For momentum- 



independent hybridization between c and / bands, this 
optical gap A opt is simply given by twice the value of the 
renormalizcd hybridization. As explained above, the hy- 
bridization is expected to scale as the square root of the 
coherence temperature, hence A opt ~ VT co hD (where D 
is the conduction-electron bandwidth) . Clearly, the sim- 
ple two-band picture falls short of capturing inelastic pro- 
cesses at non-zero energies, which will inevitably smear 
out the gap even at T — 0. Nevertheless, a pseudogap- 
like feature has been shown to survive in a(ui) when 
fully accounting for dynamic local correlation effects in 
the framework of DMFT for the standard Anderson lat- 
tice model at large U i 21 ' 27 In the results of these calcu- 
lations, the magnitude of the pseudogap has the same 
scaling as above, A opt ~ VTcoh- Remarkably, this rela- 
tion between optical gap and coherence temperature has 
been found to be nicely obeyed by a number of heavy- 
fermion metals 4^ i 10 ' 1:L i 12 However, optical conductivity 
studies in CeMIn 5 (M=Ir,Rh or Co) show little signa- 
tures of a well-defined hybridization gap.— As we show 
below, such a behavior is in principle consistent with a 
strongly momentum-dependent hybridization in the un- 
derlying Anderson lattice model. 

The finite-frequency part of the optical conductivity 
cr(u>) can be expressed through the retarded current- 
current correlation function as 

a(oj) = -n(w + irj) (23) 

UJ 

with 



n(ia;) = -|dre^(T T jt(r)j(0)) (24) 
o 

The current operator j has to calculated as the time 
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FIG. 6: (Color online) Real part of the optical conductivity 
(Ji (lu) for a) s-wave (solid), b) extended s-wave (dotted), and 
c) d x 2_ y 2-wecve hybridization (dashed). The parameters are 
as in Fig. [2ji-c, for a discussion see text. 

derivative of the polarization operator P 

j = i[W,P] (25) 

where the definition of P includes all charged particles 
(with charge g<): 

P=2 ?i R i oJo i . (26) 

As usual, approximations to the propagators and to the 
current vertex in calculating II (uj) have to be mutually 
consistent in order to respect charge conservation (ex- 
pressed by the corresponding Ward identity). 

At this point, the electrodynamics of the heavy Fermi 
liquid requires a thorough discussion. Physically, the 
/ electrons contribute to the Fermi surface and carry 
charge. While this is plausible in the Anderson model 
picture, where the charge is naturally carried by the / 
auxiliary particles, the Kondo case is more subtle: The 
/ particles of the mean-field theory are neutral spinons, 
which will carry a physical electric charge only upon in- 
clusion of gauge fluctuations, see Sec. [V] Hence, we shall 
take the Anderson model viewpoint here. In the spirit of 
the mean-field theory, we demand the current correlator 
to be calculated as the bare bubble. An expression for 
the current vertex, which is consistent with the mean- 
field propagators, is obtained from: 

jMF = i [^alm.mf, Pmf] (27) 

where charge-carrying particles Cj and fi are contained 
in P M f: 

Pmf = £ R* (4<*r + Ufa) ( 28 ) 

■ia 

Evaluating Eq. (|2T|) leads to 

Jmf - (( V k £ k)cL c k ff + ?-(V k V k )(cj cCT / kCT + &.c.)) • 

(29) 



Let us pause to emphasize that a current operator de- 
rived from Halm before the mean-field approximation 
would have a / electron contribution different from that 
in Jmf, but such a current vertex would be inconsistent 
when used together with the bubble of mean-field propa- 
gators, i.e., vertex corrections would become important. 
We point out that jmf has several shortcomings because 
Pmf treats the / as real electrons, nevertheless, the ex- 
pression (|29p is the only current operator suitable within 
the mean-field treatment of the Anderson model. We also 
note that the second term in jmf vanishes in the conven- 
tional case of a constant hybridization Vk, and ambigu- 
ities regarding the treatment of the / electrons do not 
arisej 21 ' 28 

Using the current operator jMF (|29|) . the expression for 
the real part of the optical conductivity o~\ (still assuming 
n c < 1) reads 

ai H = - V , nF{Z2k \ 2 A*5(z lk z 2k - W ) (30) 

using the abbreviations e"/ = e/ — A and 

Ak = (Vk£k) 2 (^2k - e/)(zik - ?/)+ 

+ r 2 (V k Vk) 2 ((z lk - z 2k ) + 4r 2 V k 2 ) + (31) 
+ 2r 2 V k Ve k VV k (z 2k + z lk - 2e f ) 

The result for a(uj>0) from a numerical evaluation of 
Eq. (|30)) is depicted in Fig. [6] for the three hybridization 
symmetries of Eq. (fT5|) and parameters as in Fig. [2^- 
c. In all situations, a finite gap A opt is visible in a(u>), 
which corresponds to the minimal direct gap between the 
two bands zi k and z 2 k- In the cases of the s-wave and 
extended s-wave hybridization, A opt is given by 2b and 
| Ao — lAb/^Jb 2 + t 2 , respectively; both expressions trans- 
late into A opt ~ y/TcohD (up to prefactors), as known 
before. For d ;c 2_ y 2-wave hybridization, the two bands 
cross along the nodal lines, but this crossing is at a finite 
energy away from the Fermi level. Hence, the direct gap 
is finite and given by the renormalized / level position, 
| Ao | — this translates into A opt oc T co h. Above this thresh- 
old energy, the optical conductivity follows ^Juj — A opt , 
see App. [C] 

As discussed above, a hard gap in o~{u>) will not survive 
beyond mean-field, but we expect the qualitative result 
to remain valid. We therefore conclude that a hybridiza- 
tion Vic with momentum-space nodes leads to transfer 
of optical spectral weight from the energy scale \/T co hD 
to the scale T co h (when compared to the case of con- 
stant hybridization). For actual experiments this likely 
implies that no hybridization gap will be visible in the 
optical conductivity, due to the finite width of the Drude 
peak. Such a scenario is qualitatively consistent with the 
optical-conductivity data obtained on CeMIns. 7 
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E. Thermal transport 

Low-temperature d.c. transport quantities are in prin- 
ciple candidates to probe strong anisotropies in momen- 
tum space. As an example, let us consider the thermal 
conductivity (which sometimes shows less sample depen- 
dence than the electrical conductivity). The energy cur- 
rent operator in the mean-field approximation reads^ 



Jt = ( z ik(V k 2:ik)7i T k7ik + (1 <->• 2) 



(32) 



where 7J 2 are the operators creating a quasi-particle in 
the z% t 2 band. ^From the Kubo formula, one derives the 
low-temperature thermal conductivity in relaxation-time 
approximation: 

'" = ^y] Z 2k(Vk22k)i(Vk^2k)j {-n F (z 2 k)) =T~ r 
J - — t l z 2kj 



ka 



(33) 

where T denotes the impurity-induced quasiparticle scat- 
tering rate, and we have again assumed that only the 
band z 2 k crosses the Fermi level. 

As already discussed by Moreno and Coleman 1 — in the 
context of gap-anisotropic Kondo insulators, the thermal 
conductivity will be strongly anisotropic for 3d systems 
where the hybridization has e.g. line nodes. In contrast, 
in the 2d case of a d x i_ y i hybridization, the conduc- 
tivity tensor does not have enough degrees of freedom 
to reflect the anisotropy, as the two principal axes are 
equivalent here. (The sign of /3k does not enter.) The 
same applies to hybridization function Vk with higher 
angular momenta I. Hence, for 2d anisotropic systems, 
higher-order correlation functions need to be considered, 
as e.g. probed by angle-dependent magnetoresistance; 
this is beyond the scope of this paper. [An exception 
is a p-wave hybridization (i.e. I — 1) which explicitly 
breaks the C4 rotation symmetry down to C 2 , leading to 
an in-plane transport anisotropy. Note that such a hy- 
bridization will be accompanied by a corresponding lat- 
tice distortion, which will be reflected in the entire band 
structure.] 

Finally, we note that, independent of possible trans- 
port anisotropies, the Wiedemann-Franz law will always 
be obeyed (assuming elastic scattering only) : The Lorenz 
number L, formed from the thermal conductivity Ku and 
the electrical conductivity o~n via L = k/(uT), will ap- 
proach the constant Lq = (7r 2 /3)(fcs/e) 2 in the low- 
temperature limit. This is consistent with the fact that 
we are describing a Fermi liquid. As a corollary, the re- 
cently observed violation of the Wiedemann-Franz law 
in CeCoIns at its field-induced critical points is likely 
related to inelastic scattering processes. 



V. BEYOND MEAN-FIELD THEORY 

So far, we have discussed the low-temperature proper- 
ties of "unconventional" heavy Fermi liquids using slave- 



boson mean- field theory. In principle, corrections to 
mean-field theory can be systematically taken into ac- 
count, by considering fluctuations around the saddle 
point. For the Kondo model, the correct implementa- 
tion of the Hubert space constraint, together with phase 
fluctuations of the boson field, lead to a theory where 
/ and b particles are minimally coupled to a compact 
U(l) gauge field. The Fermi-liquid phase corresponds to 
the Higgs/confining phase of the gauge theory, it is sta- 
ble w.r.t. fluctuation effects, their main effect being to 
endow the /-particle with a physical electric chargei 30 i 31 

To treat the full crossover from energies or temper- 
atures above Xk to those below T co h, different methods 
need to be employed. Local correlations can be efficiently 
captured by dynamical mean-field theory (DMFT)3 If 
DMFT is formulated for the Anderson lattice model (JTJ) , 
correlation effects arise from the local Hubbard interac- 
tion U, and consequently DMFT can be used to treat an 
Anderson model with non-local hybridization Vk as well. 
The DMFT self-consistency equation then reads: 



Galm, 100(2) — 



1 



k z-e/-E/(z)- 
I 

7 f - A(z) - V f {z) 



(34) 

= Gsiam(^) 



Here, S/ is the so-called interaction self-energy arising 
from U, and A denotes the effective hybridization func- 
tion defined by the second line of Eq. (f3~4")l . While we shall 
not numerically solve the DMFT problem (|34)) here, we 
can briefly discuss a few properties. Most importantly, 
the momentum dependence of the arising effective hy- 
bridization is dictated by the bare Vk. This implies that 
all qualitative statements in Sec. IIVI remain valid, in par- 
ticular all local moments will be fully screened at low T. 
(Technically, the DMFT reduces the lattice model (HJ to 
an effective single-impurity model, with a bath having a 
finite density of states at the Fermi level - this implies a 
fully developed Kondo effect as T — ► 0.) 

Cluster extensions of DMFT allow to handle 
momentum-dependent self-energies. Then, in principle 
the momentum dependence of the effective hybridization 
will differ from that of the bare Vk- However, we do 
not expect qualitative changes of the low-temperature 
physics described above. 

Let us note one caveat: While calculating thermo- 
dynamics and single-particle properties within DMFT 
for the Anderson model (fT]) is straightforward, electric 
transport is not. The reason is that the current opera- 
tor inevitably involves contributions from the non-local 
hybridization, see discussion in Sec. IIV Dl As a result, 
vertex corrections do not vanish in the DMFT limit, in 
contrast to standard DMFT applications^ 
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VI. TEMPERATURE-DEPENDENT 
RESISTIVITY 

In this section, we touch upon electronic properties at 
elevated temperatures. In particular, we want to focus 
on the electrical resistivity p(T) of heavy fermions with 
unconventional hybridization for T > Ik- 

In the conventional hcavy-fermion picture, the electri- 
cal resistivity p{T) at high temperatures, T 3> Ik, is 
small (ignoring phonons here), and p(T) rises upon low- 
ering the temperature due to increasing magnetic scat- 
tering. At a scale which is often identified with the lat- 
tice coherence temperature T co h, p(T) reaches a maxi- 
mum and then drops upon further cooling, behaving as 
p(T) = po + AT 2 at low T. At elevated temperatures, 
T > 7k, the scattering can be accessed using pertur- 
bation theory in the Kondo coupling, i.e., the physical 
picture is that of c electrons (with a small Fermi surface) 
scattering inelastically off the / moments. 

Bare perturbation theory gives a single-particle scat- 
tering rate 

r^aJ^l+^ln^). (35) 

A few remarks are in order: (i) In the paramagnetic phase 
of a Kondo lattice, all contributions to the conduction- 
electron self-energy up to order Jq arise from single- 
impurity scattering, (ii) The prefactor /3£ comes from 
the two external lines of the self-energy diagrams whereas 
the internal momentum summations average out all other 
form factors - this is also true for higher-order diagrams. 
Assuming that scattering arises from the local moments 
only, the simplest approximation for the conductivity 
a oc 1/p yields 

o-ij(T) cc j — -p-ui(k)u,-(k)-nk (36) 

FS 

where i>i(k) — de^/dki is the quasiparticle velocity. The 
result (|35l36p is interesting, as it shows that for form fac- 
tors /3k with nodes, the Kondo scattering is insufficient 
to render the conductivity finite, because diverges at 
least like (k—k n )~ 2 near the node at k n . The physical ori- 
gin is that conduction electrons with momenta at the hy- 
bridization nodes are not scattered at all, and this short- 
circuits all other processes, leading to infinite conductiv- 
ity. To obtain a finite conductivity, additional scattering 
needs to be considered, namely electron-electron scat- 
tering among the conduction electrons, electron-lattice 
scattering, or scattering off static impurities. The result- 
ing interplay of scattering mechanisms can be complex 
and can even modify the basic temperature dependence 
of jo(T), but we shall not analyze it here in detail. 

The physical conclusion is that the electrical current, 
at least in the temperature regime T > Tk, is primarily 
carried by conduction electrons with weak hybridization 
to the / moments, hence "nodal" quasiparticles dom- 
inate the electric transport. Recall that, in contrast, 



the low-temperature thermodynamics is dominated by 
"antinodal" quasiparticles. 



VII. COMPETITION BETWEEN KONDO 
SCREENING AND ORDERING 

As already discussed by Doniach^22 the phase diagram 
of heavy-fermion metals is determined by the competition 
between Kondo effect and inter-moment exchange (either 
of direct or RKKY type). Inter-moment exchange can 
drive magnetic ordering, but may also lead to non-trivial 
metallic spin-liquid states and to magnetically mediated 
superconductivity. The competition with Kondo screen- 
ing may be simply understood by stating that the local / 
moments can either form Kondo singlets with the conduc- 
tion electrons, or they can order in a symmetry-breaking 
fashion or else pair into inter- moment singlets. 

Thinking about these competing tendencies in momen- 
tum space, it is conceivable that, in a situation with 
momentum-space differentiation of electronic band prop- 
erties, certain ordering phenomena are favored or disfa- 
vored by a given form of the hybridization. This idea 
will be illustrated in this section, using magnetically me- 
diated superconductivity as an example, where one can 
expect an intricate interplay between hybridization and 
pairing symmetries. Concrete calculations will be per- 
formed for a Kondo-Heisenberg model in a mean-field 
approach: As in Ref. l33l a magnetic interaction between 
the / moments can be decoupled in the particle-particle 
channel, leading to pairing of spinons, which, if coexisting 
with Kondo screening, leads to BCS-type superconduc- 
tivity. 

A. Kondo-Heisenberg model 

The Anderson and Kondo lattice models (|H2p con- 
tain the competition of Kondo and RKKY interactions. 
However, in the slave-boson approach the effect of the 
RKKY interaction is lost. For mean-field calculations it 
is thus convenient to introduce an explicit inter-moment 
exchange interaction of Heisenberg type: 

Hh = ^^ s .. s ._ (37) 

The physics of the model 7Yklm + "Hh, commonly re- 
ferred to as Kondo-Heisenberg model, has been exten- 
sively discussed in the literature. We shall give a compre- 
hensive discussion of all phases and phase diagrams, but 
instead concentrate on the possible emergence of super- 
conductivity due to / electron pairing. A general frame- 
work has been laid out in Ref. 33.34 which considered 
a scenario where dominant RKKY interaction does not 
lead to antiferromagnetism, but instead to a metallic spin 
liquid state. This state, arising e.g. from geometric frus- 
tration of the inter-moment exchange, has been dubbed 
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"fractionalized Fermi liquid" (FL*), as it features light 
conduction electrons, forming a Fermi liquid, which co- 
exist with a fractionalized spin liquid formed out of the 
/ electrons. Then, in the generalized Doniach phase di- 
agram, the heavy Fermi liquid (FL) is separated from 
FL* by a quantum critical point where Kondo screening 
breaks down, but no local symmetries are hrokeni 33 i 34 
(This quantum critical point has been discussed in rela- 
tion to unconventional quantum criticality in materials 
like CeCug-xAttj; and YbRh 2 Si 2 .) 

A specific realization of FL* is a state with paired 
spinons and an emergent Z 2 gauge structure. As detailed 
in Ref. [Hi, one can expect magnetically mediated super- 
conductivity close to the quantum critical point between 
FL and a Z 2 FL*. All resulting low-temperature phases 
can be conveniently described in a mean-field approach, 
where the standard slave-boson description of the Kondo 
effect is combined with a Sp(2iV) mean- field treatment of 
the Heisenberg exchange. 35 Below, we shall extend this 
mean-field theory to the case of momentum-dependent 
hybridization. 



B. Mean- field theory and magnetically mediated 
superconductivity 

A mean-field theory for the Kondo-Heisenberg model, 
Hklm + Wh, involves a decoupling of the Kondo inter- 
action as in Sec. IIIII and of the inter-moment Heisen- 
berg exchange TLh- Using the pseudofermion represen- 
tation of the local moments as above, non-local spinon 
pairing is described by a mean field of the form Ay — 

~ (fi^fji ~ filftf)- Then, the Heisenberg interaction can 
be written at the mean-field level as 3 ^ 



H 



H.MF — 



I Ay | 2 



For time-reversal invariant states, the bond field Ay = 
Aji can be chosen to be real. Importantly, the 
mean-field Hamiltonian (f38|) is the exact solution of 
the Heisenberg model in the symplectic Sp(27V) large- 
TV limit, with a fully antisymmetric representation of 
the local moments* 3 ^ (This laige-N limit uniquely se- 
lects the particle-particle decoupling of the Heisenberg 
interaction^) Physicswise, non-zero Ay creates a para- 
magnetic phase out of the / moments; in particular, uni- 
form Ay describes a gapped Z 2 spin liquid, but also 
states with broken translational symmetry can occur 
which can be classified as valence-bond solids^ A con- 
sistent Sp(2iV) mean-field treatment of the full model 
Hklm +Hh is obtained by also decoupling the Kondo in- 
teraction in the particle-particle channel. However, one 
can show that for the Kondo part both particle-hole and 
particle-particle decoupling schemes are equivalent re- 
garding physical observables, provided that rif = 1 and 
time-reversal symmetry is present. 

The full mean-field theory is now given by Hklm.mf + 



Hh,mf, with the two "order parameters" b and A. Re- 
stricting ourselves to states without translational sym- 
metry breaking, the following mean-field phases occur: 
At high temperatures, a trivial decoupled phase with 
b = A = is realized. If the Heisenberg exchange Jh 
dominates over the Kondo coupling Jo, then A will be 
finite and b zero at low T, resulting in decoupled c and 
/ electron subsystems: This is the fractionalized Fermi 
liquid FL* described above. On the other hand, non-zero 
b and vanishing A describe a conventional heavy Fermi 
liquid FL (which was the subject of Sec. IIII[) . Finally, 
if both A and b are non-zero, Kondo screening coexists 
with spinon pairing, which leads to a true superconduct- 
ing state (SC), with pairing mediated by the magnetic 
coupling among the / moments. At sufficiently low T, 
the FL phase is always unstable towards superconduc- 
tivity in the presence of a non-zero Jh- (Note that the 
FL* phase is not a superconductor, as the / particles 
do not carry a physical charge in the absence of Kondo 
screening.) Fluctuation corrections to mean-field theory 
will smear out the finite-temperature transitions of the 
FL and FL* phases (the latter only in d = 2 dimensions), 
whereas the superconducting transition remains a true 
phase transition^ 

At this point, a more detailed discussion of the spa- 
tial structure of the Heisenberg interaction, described by 
Jn,iji and of the resulting pairing is needed. For nearest- 
neighbor exchange on the square lattice of / moments, 
each unit cell contains two bond variables, Ay . A numer- 
ical solution shows that two types of saddle points exist 
(provided that translational and time-reversal invariance 
are imposed 3 ^), namely a uniform (or "extended s-wave" ) 
solution with Ay — A on all links, and a "d I 2_ !/ 2-wave" 
solution with Ay = ±A on horizontal and vertical links, 
respectively. We will also consider the case where the 
Jn.ij only act on next-neighbor diagonal bonds, which 
together with an alternating structure of the Ay leads 
to a "cfcy-wave" mean-field solution of the Heisenberg 
part. All cases can be written in momentum space as 

Hh,mf = ^ k (ftj-kl + h - c ) + j hN& 2 (39) 



with the abbreviation Wk = — JffAak, and contains 
the "form factor" of the spinon pairing: 



cos k x + cos k y extended s-wave 
«k = ^ cos k x — cos k y d 2; 2_j / 2-wave 
2 sin k x sin k y d xy -w&ve 



(40) 



The self-consistency equation which supplement the 
Hamiltonian Hklm,mf + Hh,mf, are given by Eqs. (fH 
and 



E a k(/ kT /iki + ^-> = 2AAA. 



(41) 



All expectation values can again be expressed in terms of 
Green functions, see App. IA 21 Diagonalizing the mean- 
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field Hamiltonian, we obtain the quasiparticle energies: 

«l >2 ,3,4k = ±^= ^A 2 , + e 2 + 26 2 /? 2 + ^2 
± (-4(A e k + & 2 /£) 2 -4e£w£ 

l/2\ ^/^ 

+ (A 2 + f - k + 26 2 / 3 2 + M? 2 )^ J . (42) 

Let us point out an interesting feature of the super- 
conducting state SC: Due to the momentum-dependent 
hybridization, the internal structure of the Cooper pairs 
is highly non-trivial. In particular, the anomalous ex- 
pectation value of the physical electrons is given by 
(c k |C_ k |) oc A/?£a k M k where M k is a smooth function 
respecting the lattice symmetries (see App. IA 2|) . This 
transforms under the same representation of the space 
group as the /l k t) <x Aa k of the spinons, but has 
additional zeros from /3 k . (A somewhat similar situa- 
tion appears in the composite pairing picture of Ref . [37l . 
however, there arising from two-channel Kondo physics.) 
The additional zeros will have a strong influence on ther- 
modynamics, e.g., the power-law in the low-temperature 
specific heat will be modified. 



C. Qualitative discussion: T c <C Ik vs. T c ~ Tk 

Pairing in a Kondo-lattice system can occur in qual- 
itatively different regimes, depending on the relation 
between the superconducting T c and the characteristic 
Kondo scale. If T c <C 7k, then the proper picture is that 
of BCS-likc pairing out of a well-formed heavy Fermi liq- 
uid. In contrast, T c ~ Tk implies a strong competition 
of Kondo screening and Cooper pairing, and supercon- 
ductivity emerges out of an incoherent non-Fermi liquid 
regime. (The formal situation T c 3> Tk leads to the non- 
superconducting FL* phase.) 

Let us quickly discuss the two regimes in the framework 
of the mean-field theory, keeping in mind that inelastic 
processes at energies of order Tk will not be captured. 
The mean-field equation (|41[) reduces to an equation for 
T c if we set A to zero: 

(43) 

The factor a k originates from the fact that the momen- 
tum dependence of both the gap and the pairing interac- 
tion Wk are equal by construction. 

In the regime T c -C Tk only the quasiparticle 
band crossing the Fermi level (which we again as- 
sume to be Z2k) contributes to pairing. We replace 
the k summation by an integral over isoenergetic lines 
/ doj f u _ z c^4k/| Vz 2k |. In analogy to the standard BCS 



case we approximate tanh(u/2T c ) = 1 for ui > 2T C 
and elsewhere. The factor of a k will be taken at its 
value at the Fermi surface. Furthermore we will set 
k = (e k — z 2 k )/(z 2 k — #2k) ~ cons t: assuming that it 
is weakly varying along the Fermi surface. Neglecting 
the dependencies of the velocity and a k perpendicular to 
the Fermi surface we obtain 

A 

2 _ r du> r dA k naj ^ 
Jh J lu J 2it |V k z 2k | 

2T C FS 

with A of order the bandwidth. This gives a rough T c 
estimate of 

Tc = 2 cxp (~jW^ k K/|vz 2k |) (45) 

This equation shows the direct interplay of the form fac- 
tors: T c is enhanced if the pairing is strong (large a k ) 
in regions where quasiparticles are heavy and hence have 
large / character. In other words, "antinodal" regions 
along the Fermi surface with large hybridization are more 
susceptible to pairing. 

In contrast, in the regime of T c ~ Tk both bands con- 
tribute to pairing, and one can expect large contributions 
to the integral in Eq. (|4"3")l from the essentially flat parts 
of both bands zi,2kj which are present in particular close 
to the nodal lines of the hybridization V^. This simple ar- 
gument illustrates the competition between Kondo effect 
and pairing in the regime T c ~ Tk: "Nodal" momentum- 
space regions with less hybridization are more suscepti- 
ble to pairing - this is opposite to the statement made 
above for T c -C Tk! One should, however, keep in mind 
that the mean-field theory has limited relevance for the 
true physics at energies or temperatures of order Tk, 
the key point being that electrons in the "antinodal" re- 
gions are rather incoherent. The emerging problem of 
the pairing of incoherent fcrmions is of fundamental rel- 
evance and heavily debated for instance in the field of 
high-temperature superconductors, but rather little solid 
knowledge exists about this highly interesting strong- 
coupling phenomenon. 

D. Extended phase diagrams 

The discussion of the last section suggests that a cer- 
tain hybridization symmetry can favor or disfavor a cer- 
tain pairing symmetry. To follow up on this idea, we 
have determined mean-field phase diagrams from a fully 
self-consistent numerical solution of Eqs. (fP2|) and (|4ip . 
While our results in principle support the above state- 
ment, they also show that microscopic details of band 
structure, band filling, and pairing interaction are im- 
portant in determining T c (which may render simplistic 
arguments invalid). 

Sample phase diagrams as function of temperature T 
and Kondo coupling Jq 7 keeping Jh , t, and n c fixed, are 
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FIG. 7: Phase diagrams in dependence of temperature T 
and Kondo coupling Jo at n c = 0.3 for an extended s- and a 
d x 2_ y 2-wa,ve hybridization, combined with different types of 
a superconducting symmetry: a) d x 2_ y 2-w&ve, b) extended 
s-wave, and c) d xy -wa,ve. The left-panel insets show the real- 
space structure of spinon pairing fields A„, leading to the 
certain type of pairing symmetry. Thick (thin) lines refer to 
first (second) order phase transitions, for further details see 
text. 



shown in Fig. for hybridizations of extended s- and 
d x 2_ y 2-tyi>e and various pairing symmetries. The overall 
structure of the phase diagram was discussed above in 
Sec. I VII Bl and is identical to that described in Ref. l33l . 

The different pairing symmetries in Fig. [7J have to be 
understood as follows: For a given Heisenberg interac- 
tion, saddle points with different spinon pairing sym- 
metry occur, and to plot the phase diagrams we have 
restricted our attention by hand to one of the saddle 




-0.955 



-0.96 




FIG. 8: (Color online) Free energy plotted against the Kondo 
coupling Jo at n c = 0.3 and T = 10 -5 for extended s-wave hy- 
bridization (left) and d x 2_ y 2-wave hybridization (right) and 
different superconducting symmetries: a) d x 2_ y 2-wa,ve, b) ex- 
tended s-wave, and c) d xy -wa.ve (compare to the correspond- 
ing phase diagrams shown in Fig. [7J The thin dashed line 
shows the normal-state solution with A = 0. 



points. The correct superconducting phase is obtained 
by comparing the free energies, given by the mean-field 
expression 



E(^-Ao) 



(46) 



- A/>n c + 7VA + MJ H A 2 + N— . 

Jo 



Plots of the free energies at low T are presented in Fig. 

m 

The rough conclusion for the particular dispersion and 
band filling used here is that extended s-wave hybridiza- 
tion favors the d x 2_„2-wave superconducting symmetry, 
while d a; 2_j / 2-wave hybridization favors d xy -wave super- 
conductivity. Given the structure of the Fermi surfaces 
in the FL phase, this is not unexpected: If both and 
/3k have extended s structure, the pairing interaction is 
rather small near the Fermi surface (Fig. \2)p) , and con- 
sequently the superconductivity is weak in the left panel 
of Fig. [7Jd, whereas both <i-wave pairing states perform 
nicely in energy. For /3k of d x 2_ y 2 form, the Fermi surface 
(Fig. [2b) is mainly located close to the momentum-space 
diagonals, favoring d xy pairing (right panel of Fig. [7b). 



VIII. CONCLUSIONS 

In this paper, we have explored the consequences 
of a strongly momentum-dependent hybridization be- 
tween conduction and local-moment electrons in heavy- 
fcrmion metals: In the Fermi-liquid regime, the quasi- 
particle properties become strongly anisotropic along the 
Fermi surface: "nodal" quasiparticles are light c elec- 
trons, whereas "antinodal" quasiparticles are heavy and 
have essentially / character. An interesting dichotomy 
arises: While the low-temperature thermodynamics is 
dominated by heavy antinodal quasiparticles, the elec- 
trical conductivity at elevated temperatures is carried by 
unhybridized nodal quasiparticles. Experimentally im- 
portant is the low-temperature optical conductivity c(oj): 
Due to the strongly momentum-dependent gap between 
the effective bands, the hybridization gap in u{lo) is es- 
sentially smeared out. 

Further, we have advocated the idea that the 
momentum-space structure of the hybridization is impor- 
tant in selecting ordering phenomena which compete with 
Kondo screening near quantum criticality. Here, two 
regimes need to be distinguished: For energies or temper- 
atures T much smaller than the coherence temperature 
Tcoh, a weak-coupling quasiparticle picture is often ap- 
propriate, and instabilities of the heavy Fermi liquid are 
determined by the interaction among the (anisotropic!) 
quasiparticles. In contrast, for T ~ T co h fascinating 
strong-coupling phenomena can be expected, for example 
unconventional superconductivity emerging from a non- 
Fermi liquid regime. This physics will be dominated by 
inelastic processes, which again are strongly anisotropic 
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in momentum space. A detailed study should be under- 
taken using cluster extensions of dynamical mean-field 
theory, but is beyond the scope of this paper. 

On the experimental side, CeNiSn and CeRhSb have 
been established to be half-filled Kondo semimetals with 
a hybridization gap vanishing along a certain crys- 
tallographic axis . 13 ' 14 i 39 The CeMIns compounds are 
candidates for Kondo metals with strongly anisotropic 
hybridization^ but other Ce or Yb materials where a 
clear-cut hybridization gap in ct(oj) is absent may fall 
into this category as well. We note that first-principles 
calculations based on density- functional theory could, in 
principle, be able to determine the hybridization symme- 
try, but strong interaction effects can render the conclu- 
sions invalid. Recent x-ray absorption studies are promis- 
ing in paving a way to an experimental determination 
of the required microscopic information^ To probe the 
anisotropic quasiparticle properties in the Kondo regime, 
high-resolution angle-resolved photoemission is the ideal 
tool (with the restriction that in can only be applied 
to quasi-2d systems). As outlined in Sec. IVI1 unusual 
behavior in the finite-temperature resistivity may also 
be connected to nodes in the hybridization function. 
Clearly, more detailed theoretical investigations of trans- 
port properties are needed. Finally, we mention that the 
strong-coupling pairing regime T c ~ Tk is very likely re- 
alized in the fascinating superconductor PuCoGas.— 
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Notes added. While this paper was being completed, 
a related paper by Ghaemi et al. i2 on angle-dependent 
quasiparticle weights appeared on the arXiv. Their re- 
sults for the low-temperature regime of anisotropic heavy 
Fermi liquids are related to ours. 

After submission of this paper, Shim et al. 43 published 
a paper on first-principles calculations for Celrlns, which 
support the idea of a strongly momentum-dependent hy- 
bridization function. 



1. Green functions 

The Kondo-lattice mean-field Hamiltonian can be 
rewritten in a matrix form: 



t ( £k b(3 k 
bfa -A 



<IV + const. 



(Al) 



with vffjc = ^c kcr , . In the following, we shall denote 

retarded Green functions 



Gab(z) 



die 1 



(-ie(t){[A(t),B(0)\. 



(A2) 



as ((A:B)) Z . Defining the matrix propagator 

we obtain by explicit inversion 
G(k,z) 



(A3) 



_ ( (( c L;0)* (( c L;/k,T» z 



((/kcr' C kcr)) z ((/kcr! /kcr)) ; 

1 fz + Xo bp k 

(z - Zl k )(z - Z 2 k) V Z " ik 



(A4) 



The thermal expectation values required for the mean- 
field equation are obtained by summing over Matsubara 
frequencies; this can be done analytically, as the excita- 
tion energies, Eq. (fl7|) , are known. 



2. Green functions in the presence of a Heisenberg 

term 

The Hamiltonian containing the additional Heisenberg 
term, Eq. (|39|) . has to be rewritten in a matrix form in 
analogy to App. IA 11 The inversion of (z — H^) pro- 
vides the needed Green functions. We use the shorthand 
h{z) = Y[ l (z-z i ). 



((/kT' /k|)) z — 



(z - e k )((z - \ ){z + e k ) - b 2 j 2 ) 



Hz) 



(A5a) 



7\ _(z + e k )((z + \ Q )(z - e k ) - b 2 -y 2 ) 



((/-kj.; /-kj.))z ~ 



Hz) 



(A5b) 



APPENDIX A: MEAN-FIELD THEORY 

In this appendix, we list the expressions of Green func- 
tions required for the implementation of the mean-field 
theory. 



((/kT i /l k [ )) 2 — ((f-kl'fk'\)) z ~ 



kj.' Jky 



H z ) 



«4;<£ k ,»« = «c- k ,;0>. = (A5d) 



Hz) 
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<<4 T ;/ k T>>*H</^ c kT>> 



We can see that the mean-field equations of both theories 
b/3k((z — \q)(z + 6k) — b 2 /3 k ) (A5e) correspond to each other, and rV and b are the effective 
JJ^) band hybridizations in the two-band model. 



«/- ki ;c t - ki )) z = « C _ ki ;/i ki )) 2 

_ -6/3 k ((z + A )(z-e k )-6 2 /3 k ) 



h(z) 



(A5f) 



( ( c kt ' C k|)) z — 



(( c - k i; c - k i))^ 



-fr 2 /? 2 (z + Ap) + (z + e k )(z 2 - VF 2 - Aq) 

(A5g) 

_ -& 2 /3 2 (z - A ) + (* - e k )(z 2 - ^ k 2 - Aq) 



h(z) 



(A5h) 



APPENDIX B: EQUIVALENCE OF ANDERSON 
AND KONDO MEAN-FIELD THEORIES 

Here we compare the two sets of mean-held equations 
for the Anderson and Kondo lattice models, Eqs. © 
and (fl"2|) . These are expected to be equivalent, once the 
Kondo limit is taken in the Anderson-model equations. 
For U — ¥ oo, the Kondo coupling is Jo = V 2 /|e/|. Fur- 
ther, the Kondo limit implies r — > 0, for otherwise the 
effective hybridization would diverge. The average / oc- 
cupation then becomes unity, and the / and / operators 
are equivalent. The physical valence fluctuations in the 
Anderson model are projected out by e/ — > — oo. In this 
limit, the effective / level energy, (e/ — A), stays finite (i.e. 
a fraction of the bandwidth) to ensure Y^aifafv) = 1- 
Therefore, V 2 /\X\ -> J Q . 

With this knowledge about the limiting behaviors, the 
first of the Anderson-model mean-field equations (JHJa,) 
transforms like 



rV 



V 2 (3 k 



2Af 



/ k <T C k <T + C kafka )) 



ka 



(Bl) 



APPENDIX C: OPTICAL CONDUCTIVITY 

We briefly discuss the behavior of the inter-band op- 
tical conductivity close to the threshold energy, for the 
case of d- wave hybridization. In Eq. (|30p , the matrix ele- 
ments are non-singular near the threshold, hence we only 
need to analyze the behavior of 



V /(A + e k ) 2 +46 2 /3 I 



(CI) 



A change of variables is a: 



fox */2 x 



2k ) 



(C2a) 



V2 



( k x + ky) 



(C2b) 



fco is determined by e(fco,fco) = For u> near — Ao, 
AE(k x , k y ) can be approximated by a second-order Tay- 
lor expansion in k x and k y around 0. The k integra- 
tion is restricted to the Fermi sea because of the factor 
nF{zik)- The Fermi surface can approximated by also ex- 
panding around k x , k y = 0. Power counting in the inte- 
gral then shows, that the first non- vanishing contribution 
to the optical conductivity for a d-wave hybridization is 

CX \JuJ + Aq. 
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